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Quantum state filtering is a variant of the unambiguous state discrimination problem: the states 
are grouped in sets and we want to determine to which particular set a given input state belongs. The 
simplest case, when the N given states are divided into two subsets and the first set consists of one 
state only while the second consists of all of the remaining states, is termed quantum state filtering. 
We derived previously the optimal strategy for the case of N non-orthogonal states, {|^i}, . . . , |V"Jv)}> 
for distinguishing from the set {Jf/^), . . . , |i/>jv)} and the corresponding optimal success and 
failure probabilities. In a previous paper [Phys. Rev. Lett. 90, 257901 (2003)], we sketched an 
appplication of the results to probabilistic quantum algorithms. Here we fill the gaps and give the 
complete derivation of the probabilstic quantum algorithm that can optimally distinguish between 
two classes of Boolean functions, that of the balanced functions and that of the biased functions. 
The algorithm is probabilistic, it fails sometimes but when it does it lets us know that it did. Our 
approach can be considered as a generalization of the Deutsch-Jozsa algorithm that was developed 
for the discrimination of balanced and constant Boolean functions. 

PACS numbers: 03.67-a 



I. INTRODUCTION 

The Deutsch-Jozsa algorithm was one of the first quan- 
tum algorithms 1] . It makes it possible to perform on a 
quantum computer in one step a task that, on a classi- 
cal computer, would require many steps. In particular, 
one is given a "black box" that, when given an n bit 
input, returns either or 1. The box simply evaluates 
a Boolean function, that is a function that maps n-bit 
binary numbers to the set {0, 1}. It is guaranteed that 
the functions are of one of two types, they are either con- 
stant, giving the same value for all inputs, or balanced, 
giving for half the inputs and 1 for the other half. The 
task is to determine whether the function we have been 
given is constant or balanced. On a classical computer, 
we would, in the worst case, have to evaluate the func- 
tion 2(" _1 ) + 1 times to determine this. On a quantum 
computer the determination can be made with certainty 
with only one evaluation of the function. 

The quantum algorithm works by mapping the func- 
tion onto a quantum state. The quantum states for con- 
stant and balanced functions lie in orthogonal subspaces, 
and they can, therefore, be distinguished perfectly. An 
obvious question to ask is whether the method can be 
extended to distinguish functions whose vectors do not 
lie in orthogonal subspaces. Here we shall show that this 
is indeed the case, and to do so we shall make use of a 
procedure known as quantum state filtering 

Quantum state filtering considers the following prob- 
lem. We are given a system that is in one of a known set 
of N quantum states {ipi, . . . iPn}- Our task is to deter- 
mine whether the system is in ipi or not, that is, whether 
it is in the set {ipi} or the set {V>2, ■ ■ ■ iPn}- This is what 
quantum state filtering accomplishes. It is a probabilis- 
tic procedure, it may fail, but we know when it does. 
The failure probability should be as small as possible, 
and how this is accomplished depends on the a priori 



probabilities with which the states {V'l, ■ • • V'jv} appear 
and the overlaps of ipi with the other states. In order 
to generalize the Deutsch-Jozsa algorithm, we can apply 
quantum state filtering to the quantum states into which 
the Boolean functions are mapped. 

Quantum state filtering is an instance of mixed state 
discrimination, in particular the discrimination of a pure 
state from a mixed state. The set of states {V>2, ■ ■ ■ Vw} 
along with their a priori probabilities can be viewed as 
an ensemble, which can be represented by a density ma- 
trix. There are several different strategies one can adopt 
in discriminating mixed states. One can minimize the 
probability of making an error 0, [f| Q • A second strat- 
egy is that of unambiguous discrimination in which one 
never makes a mistake in identifying the state, but one 
can sometimes fail to obtain any information about the 
state one was given Quantum state filtering is an 

example of this approach. There are also hybrid strate- 
gies in which one can make mistakes and one can also 
fail to obtain an answer [l4l fl5| . The addition of the 
option of failing to obtain an answer makes it possible to 
obtain a smaller error probability than is possible when 
the possibility of failing is not present. 

The paper is organized as follows. In the next section 
we briefly review quantum state filtering, following £| 
but with a notation simplified for the present purposes. 
In Section 3 we apply state filtering to the problem of 
distinguishing balanced functions from particular biased 
functions, i.e. functions that have a preponderance of ze- 
roes or ones. In Section 4 we provide a specific example 
of the POVM that figures in the state filtering procedure, 
and we summarize our results in the Conclusion. 



II. QUANTUM STATE FILTERING 

As was mentioned in the Introduction, what we wish to 
do, is to determine whether a system we have been given 
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is in the state ip\ or not, given that it must be in one of 
the states {■0i, ■ ■ ■ iI>n} an d that the state ipj occurs with 
probability r\j . This can be accomplished in one of three 
ways, and which way is the best depends on the a priori 
probabilities and the overlaps of the states. Two of these 
methods are von Neumann projective measurements, and 
the third is a POVM. 

The first projective measurement projects onto the 
subspace orthogonal to ipi. This is accomplished by the 
projection operator 

= 1-1^x^1. (2.i) 

If we obtain the value 1, then the procedure has suc- 
ceeded, and we know that the system we have was not 
in the state ipi- If we obtain 0, then the procedure has 
failed. This happens with a probability 

Qi=rji+S, (2.2) 

where 

N 

S = 5>|<Vi|^->| 2 (2-3) 

is the average overlap between the two subsets. 

A second possibility is to split into two compo- 
nents, \ip\) — \ipi) + |^|)- Here \ipi) is orthogonal 
to the subspace, H2, that is spanned by the vectors 
1^2)) • • • \iPn), and \ipf) lies in H.2- Their normalized ver- 

sions are \ft) = \ipf}/W£\\ and |$|) = |V>!>/||V!ll, re- 
spectively. We now introduce the operator 

= - (I- |$>$!|), (2.4) 

which has eigenvalues 1,0, and —1. If we measure 
and obtain 1, then the vector was l^i)) if we obtain — 1, 
then the vector was in the set {IV^), ■ ■ ■ IV'iv)}, and if we 
obtain 0, the procedure failed. In this case the probability 
of failure, Q2, is given by 

Q 2 = ??1 ||^|| 2 + ^-. (2.5) 

Will 

Which of these two particular strategies is better is de- 
termined by which of these two failure probabilities is 

smaller. In particular, Qi > Q2 if ?7i||V'il| 2 > S, and vice 
versa. 

The third measurement is based on a positive-operator 
valued measure (POVM, 0), and it can do better in an 
intermediate range of parameters. The POVM can be 
implemented by a unitary evolution on a larger space 
and a selective measurement. The larger space consists 
of two orthogonal subspaccs, the original system space 
and a failure space. The idea is that the unitary evolu- 
tion transforms the input sets into orthogonal sets in the 
original system space and maps them onto the same vec- 
tor in the failure space. A click in the detector measuring 



along this vector corresponds to failure of the procedure, 
since all inputs are mapped onto the same output. A no- 
click corresponds to success since now the non-orthogonal 
input sets are transformed into orthogonal output sets in 
the system space. The one-dimensionality of the failure 
space follows from the requirement that the filtering is 
optimum, as shown by the following simple considera- 
tions. Suppose that ipi is mapped onto some vector in 
the failure space and the inputs from the other set are 
mapped onto vectors that have components perpendicu- 
lar to this vector. Then a single von Neumann measure- 
ment along the orthogonal direction could identify the 
input as being from the second set, i.e. further filtering 
would be possible, lowering the failure probability and, 
contrary to our assumption, our original filtering could 
not have been optimal. 

In particular, let TLs be the D-dimensional system 
space spanned by the vectors . . . \iPn)} where, ob- 

viously, D < N. We now embed this space in a space of 
D + l dimensions, TLs+A = 7~Is®7~La, where H.a is a one- 
dimensional auxiliary Hilbert space, the failure space or 
ancilla. The basis in this space is denoted by \4>a), where 
\\4>a\\ = 1- Thus, the unitary evolution on Hs+A is spec- 
ified by the requirement that for any of the input states 
(J = 1, . . . , N) the final state has the structure 

\^) out = U\1>j) = y/pjlrt) + v^e^l^), (2.6) 

where \ipj) £ Ti-s, and \\ipj\\ = 1 From unitarity the rela- 
tion, pj + qj = 1, follows. Furthermore, pj is the proba- 
bility that the transformation \tpj) — ► |?//) succeeds and 
qj is the probability that \ipj) is mapped onto the state 
\4>a)- ln order to identify pj and qj with the state-specific 
success and failure probability for quantum filtering we 
have to require that 

WK-)=0, (2.7) 

for j = 2, . . . N. We can now set up D + 1 detectors in 
the following way. One of them is directed along 
D — 1 along the remaining D — 1 orthogonal directions 
in the original system space and the last one along \4>a) 
in the failure space. If any of the system-space detectors 
clicks we can uniquely assign the state we were given to 
one or the other subset and a click in the failure-space 
detector indicates that filtering has failed. It should be 
noted that, in general, for the unitary operator given in 
Eq. (ESJ) to exist we must have D = N 

In order to optimize the POVM, we have to determine 
those values of qj in Eq. 1)2. 6fl that yield the smallest 
average failure probability Q, where 

N 

Q !>..'/,- (2-8) 

Taking the scalar product of U\tpi) and U\ipj) in Eq. 
(|2.6|) . and using Eq. (|2.7|) . gives 

\(H^j)\ 2 =H^, (2-9) 
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for j = 2, ... , N, and Eq. I|2.8[l can be cast in the form 
QiQi) — + S/<li- Unitarity of the transformation 
U delivers the necessary condition that q\ must lie in 
the range ||'0il| 2 < ?i < 1- Details of the derivation, 
along with a discussion of the sufficient conditions for 
the existence of U , can be found in Q. Provided that a 
POVM-solution exists, the minimum of Q(q\) is reached 
for <?! = y/S/r)\ and is given by 

Qpovm = 2^S. (2.10) 

Thus, the failure probability for optimal unambiguous 
quantum state filtering can be summarized as 

f 2VmS tfmU!\\ 4 <s<m, 

Q=\ Vi+S if S > rji, (2.H) 

1*11* + ifS <77l ||^|| 4 . 

The first line represents the POVM result, Eq. fFTty . 
and it gives a smaller failure probability, in its range of 
validity, than the von Neumann measurements, Eqs. (|2.2() 
and Eq. I|2.5[> . Outside of the POVM range of validity we 
recover the von Neumann results. It should be noted that 
for these results to hold, unlike for unambiguous state 
discrimination, linear independence of all states is not 
required. Instead, the less stringent requirement of the 
linear independence of the sets is sufficient, in agreement 
with the findings in |l7j . 

It is useful to see what happens if the POVM is applied 
to a vector that lies in the subspace, Ti.2, spanned by the 
vectors {tp2, ■ ■ ■ i[>n}, but is not one of these vectors. Let 

JV 

M = I>hfc>> ( 2 - 12 ) 

J'=2 

so that 

UW) = Vp^+VQ^IM, (2.13) 
where \\ip'\\ = 1, and 

N 

3=2 

N 

We note that p is the probability that we would determine 
that the state ip is in Ti.2, and q is the probability that 
we would fail to do so. Taking the inner product of the 
above equation with U\ipi) gives us that 

q 1 q = \(i h \i>)\ 2 . (2.15) 

Therefore, the POVM can be used to solve a more general 
problem than the one for which it was designed. In par- 
ticular, suppose we are given a system that is guaranteed 



to be either in the state ipi or in a state in the subspace 
Ti.2, and we want to determine which is the case. We can 
do this by applying the POVM presented here, and our 
probability of failing is given by the above equation. Be- 
cause the POVM was designed for a particular basis of 
Ti.2, it will not necessarily be optimal for the more general 
problem, but it will work. 



III. APPLICATION TO A BASIS FOR 
BOOLEAN FUNCTIONS 

We can now apply this result to distinguishing between 
sets of Boolean functions. Let f(x), where < x < 2 n — 1, 
be a Boolean function, i.e. f(x) is either or 1. One of the 
sets we want to consider is the set of balanced functions. 
The other will consist of two "biased" functions. We shall 
call a function biased if it is not balanced or constant, i.e. 
if it returns on mo of its arguments, 1 on mi = 2™ — mo, 
and mo 7^ mi ^ or 2™ _1 . In order to discriminate a 
particular biased function from an unknown balanced one 
2( n_1 ) -\-mi + 1 function evaluations are necessary in the 
worst case classically, where we have assumed, without 
loss of generality, that mi < mo. This number comes 
from the fact that there are balanced functions that agree 
with our biased function in 2(™- 1 ) +mi places. Suppose 
that a balanced function is 1 for every argument for which 
our biased function is 1. This means that they agree for 
at least mi arguments. Of the remaining arguements, the 
biased function will be on all of them, and the balanced 
function will be on 2(™~ 1 ) of them. This means that 
the functions agree in a total of 2 ( ™~ 1 ) + mi places. 

As has been mentioned, the second set of functions we 
shall consider has only two members, and we shall call it 
Wk. A function is in W k if f(x) = for < x < [(2 fe - 
l)/2 fe ]2" and f(x) = 1 for [{2 k ~ l)/2 fc ]2™ < x < 2 n - 1, 
or if f(x) = 1 for < x < [{2 k - l)/2 fc ]2" and f(x) = 
for [(2 k - l)/2 fc ]2" < x < 2 n - 1. The problem we wish 
to consider is distinguishing between balanced functions 
and functions in Wk, that is, we are given an unknown 
function that is in one of the two sets, and we want to find 
out which set it is in. We note that the two functions in 
Wit are biased functions, so that this problem is a partic- 
ular instance of a more general problem of distinguishing 
a particular set of biased functions from balanced func- 
tions. This is by no means the only example the method 
we are proposing here can handle, but it is a particularly 
simple one. 

This is clearly a variant of the Deutsch- Jozsa problem 
1] . In that case one is given an unknown function that is 
either balanced or constant, and one wants to determine 
which. Classically, in the worst case one would have to 
evaluate the function D/2 + 1 times, where we have set 
D = 2", but in the quantum case only one function eval- 
uation is necessary. The solution of this problem makes 
use of the unitary mapping 

\x)\y)^\x)\y + f(x)), (3.1) 
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where the first state, |a;), is an n-qubit state, the second 
state, \y), is a single qubit state, and the addition is 
modulo 2. The state |x), where x is an n-digit binary 
number, is a member of the computational basis for n 
qubits, and the state \y), where y is either or f, is a 
member of the computational basis for a single qubit. 
In solving the Deutsch-Jozsa problem, this mapping is 
employed in the following way 



1 



D-l 



£ l*>(|o>-|i» 



i 



D-l 



£(-l)/C*)|z)(|0>-|l>), 



V2D V2D 

(3.2) 

and we shall do the same. This has the effect of map- 
ping Boolean functions to vectors in the D-dimensional 
Hilbert space, Hd\ the final qubit is not entangled with 
the remaining n qubits and can be discarded. The vectors 
Ylx=o(~ iy^\%) that are produced by balanced func- 
tions are orthogonal to those produced by constant func- 
tions. This is why the Deutsch-Jozsa problem is easy 
to solve quantum mechanically. In our case, the vectors 
produced by functions in Wk are not orthogonal to those 
produced by balanced functions. However, we can solve 
the problem probabilistically by using our unambiguous 
state filtering procedure. 

In order to apply quantum state filtering to this prob- 
lem, we note that both functions in Wk are mapped, up to 
an overall sign, to the same vector in Hd , which we shall 
call \wk). The vectors that correspond to balanced func- 
tions are contained in the subspace, Hb, of Hd, where 

Hb = {\v) E Hd\Y1x=o v % = °}> and v * = ( x \ v )- Tnis 
subspace has dimension 2™ — 1, and it is possible to choose 
an orthonormal basis, = 2, . . . D}, for it in which 

each basis element corresponds to a particular balanced 
Boolean function. As any state in Tib is a linear combina- 
tion of the \vj) basis vectors, a state filtering procedure 
that can optimally distinguish \wk) from the set 
will also discriminate \wk) from any set of states in Hb, 
though not necessarily optimally. 

Let us now construct the basis for Hd that we wish to 
use. This will be done inductively. For D = 2, choose 
the basis vectors to be 

\v{ 1) ) = 4(1,1) |,;«) = i=(l,-l)- (3.3) 



V2" 



Note that |i4^) is the vector that corresponds to the 
two balanced one-bit Boolean functions, and corre- 
sponds to the two constant functions. These vectors are 
orthonormal. For D = 2 2 the basis vectors are 



V2" 



l^ 2) > = 1,1,1) K0 = |(i,i,-i,-i) 

\4 2) ) = ^(1, -1, 1,-1) l4 2) > = |(1, -1, -1, !)• (3-4) 

(2) 

The vector \v\ ) corresponds to the two constant func- 
tions, and the other basis elements correspond to bal- 
anced functions. 



1 



,(2)\ - I, 



Now let us construct a basis for arbitrary D = 2™. 
First, suppose that we have all of the bases {v™\j — 
1, . . . 2 m } for m < n — 1. We now want to construct the 
basis for D — 2". Denote a string of p ones by u(p). Our 
first basis vector is just 



» 



1 



«£>)), 



(3.5) 



which just corresponds to a constant function on n bits, 
and our second is 



(u(D/2),-u(D/2)), 



(3.6) 



The general basis vector is composed of blocks of ones 
and minus ones. It can be expressed as 



(v 



(p-i). 



2M ,(«(J}/f)-«(I)/f)) 



(3.7) 



where (vj k 1 ') q is just the qih component of Vj 1 . Here 
we have that j runs from 1 to 2 P_1 , and p goes from 2 
to n — 1 . That these vectors are orthogonal for different 
values of p can be seen rather easily. Let the first vector 
have p = pi and the second have p — p2, where p\ < P2- 
When taking the inner product, each block of the form 
(u(D/2P 1 ) 1 -u{D/2P 1 ), in the first vector will be paired 
with a block of uniform ones or minus ones in the sec- 
ond vector. This means that each (u(D/2 Pl ), -u(D/2 Pl ) 
block in the first vector will contribute a zero to the total 
inner product, so that the entire inner product is itself 
just zero. When taking the inner product of two vectors 
for which the values of p are the same, we see that the in- 



(p-i)l,(p-i)\ 



where 



ner product is just proportional to (u^ 
ji and j2 are the j values for the two vectors. Therefore, 
if ji 32 the vectors are orthogonal. Henceforth, we 
shall denote the vector in Eq. 1)3. 7JI by v p j where j runs 
from 1 to 2 P_1 , and p goes from 2 to n — 1. The vector 
in Eq. (|3.5|l will be denoted as t>o,i, and the vector in Eq. 
(|3.6() as vi^. 

Let us first see how the filtering procedure performs 
when applied to the problem of distinguishing \wk) from 
the D — 1 orthonormal basis states, \v p j), for p > 0, 
in Hb- We assume that their a priori probabilities are 
equal, and we shall denote this probability by 77, where 
r/ = (1 — 7/1)/ (2™ — 1) and 771 is the a priori probability 
for the function to be in Wk- Together with the state 
\wk), the total number of states is D, the dimension of 
the system Hilbert space. 

In order to find which of the three measurement proce- 
dures is optimal, we need to calculate both S and 
Now \wk) is a unit vector, so that the sum of the square of 

its component along \vq,i) 
we have that 



and ||7i>j!|| 2 is one. Therefore, 



ul\\ 2 = 1-IK1K) 

2 k - 1 



2(2fc-2) 



(3.8) 
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The calculation of S is particularly simple in this case. 
We have 



n-1 2"- 1 

s = ^Y \( v pd\ w k)\ 2 

p=l j=l 
= r,(l - \(v 0< x\w k )\ 2 ) = 
= ■ 



T](2 k - 1) 
2(2fe"2) 



(3.9) 



Substituting these quantities into the conditions in Eq. 
(|2.11() , wc find that the first von Neumann measurement 
is optimal if r\\ > £i, the POVM is optimal if £2 < r\\ < 
£1, and the second von Neumann measurement is optimal 
if r\\ < (2 where 



Ci = 

C2 = 



(2™ - l)(2 fc - 1) 



2 2(fc-l) 

2 k - 1 



-(n-fc+2) 



2 2(fe-l)( 2 n _ 1) + 2 k 



1 



2 -(n+fe-2) ) ( 3 . 10 ) 



where, in the approximate expressions, we have assumed 
that 2™ > 2 k > 1. 

It is useful to get an idea of how much the failure prob- 
abilities of the different prodecures differ. To do so, we 
will look in the range in which the POVM is optimal, 
and compare the failure probabilities of the three differ- 
ent measurements. In this range, the failure probabilities 
are given by 



Qpovm 
Qi 
Q2 



1 



2 k-2 

Vi + 



m (l - m )(2 k - 1) 
2" - 1 

(l-^)(2 fc -l) 



1/2 



2 2k - 2 (2 n - 1) 

T?l(2 fc -1) 1-Tft 
2 2fe-2 + 2 n - 1 ' 



It should be noted that while Qpovm in the above ex- 
pression is valid only if Q2 < ??i < Ci ! t ne expressions for 
Qi and Q2 are valid for any value of 771. Assuming that 
2 n > 2 k > 1 we find that 



Qpovm 
Qi 



Qpovm 



n-\-k 



(3.12) 

so that the POVM result represents a considerable im- 
provement over either of the von Neumann measurements 
when it is optimal. For example, in the case in which all 
of the a priori probabilities are equal, i.e. r\\ = 1/2™, we 
have that 



Q 



POVM — 



{2 k -l) 1 ' 2 

2 n+fc-2 



1 



2 n-2+(fc/2) ' 



(3.13) 



both of the ratios in Eq. (|3~T2l arc 4/2 fc/2 . For k > 1, 
this implies that the difference in performance between 
the POVM and the von Neumann measurements can be 
significant. 



IV. EXAMPLE OF POVM 

So far, our discussion of the POVM has been rather 
abstract; we know when it exists and when it does not, 
and we have described how it works. It is useful, how- 
ever, to see explicitly how the unitary operator (see Eq. 
H2.6|l ) that plays a prominent role in the scheme can be 
constructed. Our task is to find explicit expressions for 
the vectors \ip'j}- In order to do this, we shall consider an 
explicit example, the case k = 2. 

We first note that (w2\v p j) = for all p > 2, which 
means that all of these vectors can be perfectly distin- 
guished from \w2}- Therefore, we define U\v Pt j) = \v p j) 
for p > 2. U will map the remaining four-dimensional 
subspace, spanned by the vectors with p < 2, into 

itself, and we can henceforth confine our attention to this 
subspace. 

Let us now set 



l^i) = \w 2 ) \ip2) = Ka) 

^3) = |«2,i) \i>) = \v2.2) 



(4.1) 



The condition that guarantees the existence of U is that 
the matrix, M , given by 



M jk - (^M> = WiM - Jw^ k 



(4.2) 



be positive. Taking the inner product of U\^>j) and 
^1^1)1 j 1) we nnd that 



so that we can express Mjk, for j, k > 1, as 

Mjk = {ipj\ipk) ~ -(^#1)^11^) 
(3.11) Setting x = l/(4gi), we find that 



(4.3) 



(4.4) 



M 



( 1 - qi 

1 — x —x x 

—x 1 — x x 

\ x x 1 — x 



(4.5) 



The eigenvalues of this matrix are 1 — (71, 1 — [3/(4q>i)], 
and 1, which is doubly degenerate. The matrix is positive 
if (3/4) < qi < 1. 

We can find the vectors \ip'j) by noting that, because 
M is positive, it can be expressed as M = B^B, where 
B = Uqt/M, and Uq is an arbitrary unitary operator. If 
we set 

|V>j-> =Bl>0,l>+ B 2>l,l)+£ 3 >2,l)+£4>2,2), (4-6) 

then we find that Mjk = (ip'j \ip'k), and we have completely 
specified U. Note that there is considerable arbitrariness, 
because we are free to choose the unitary operator Uq- A 
particular choice yields 



K) = V 1 - <?iKi), 
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1^2,2) 



ma) 



1 - 3K2) 



(4.7) 



These vectors, plus the choice of q\, which one finds by 
minimizing the failure probability and checking to see 
whether it is in the allowed range (in this case between 
3/4 and 1), completely specify the POVM. 



V. APPLICATION TO ALL BALANCED 
FUNCTIONS 

Now that we know how this procedure performs on the 
basis vectors in Tib, we shall examine its performance on 
any even function. What we shall do is compare the 
performance of the two von Neumann measurements to 
that of the POVM that is optimal for the basis vectors. 
We cannot apply the filtering procedure directly to the 
set of balanced functions, unless we want the dimension 
of the space in which the POVM is defined to be equal 
to the number of balanced functions. We are instead 
interested in a POVM acting in the space Tin © 'Ha, 
where TIa is a one-dimensional auxiliary Hilbert space. 
This can be accomplished by using the POVM that was 
derived for the more restricted problem of distinguishing 
the basis vectors from \iVk)- 

First let us see how the von Neumann measurements 
perform. We shall assume that \wk) has an a priori prob- 
ability of rji and that each of the vectors corresponding 
to a balanced function has the same probability. There 
are 



bat 



D 
D/2 



(5.1) 



vectors corresponding to balanced functions, so that the 
probability of each of them is (1 — r\\)jMb a i- In order 
to calculate the average error probabilities we need to 
calculate 



S b 



M ba l 



(5.2) 



where the sum is over all vectors corresponding to bal- 
anced functions. This sum is calculated in the Appendix, 
and we find that 



S h = 



1 - Vi 
D-l 



(5.3) 



where = (2 fc — l)/2 2fc ~ 2 . This expression is identical 
to the one calculated for the basis vectors alone. That 
means that the failure probabilities for the two von Neu- 
mann measurements, in the case of all balanced functions 



and not just the basis vectors, are still given by the ex- 
pressions in Eq. (|3.11l) . 

Now let us turn our attention to the POVM. As dis- 
cussed at the end of Section 2, even for the POVM de- 
signed to distinguish the basis vectors from \ wk), the fail- 
ure probability for any balanced function vector, \vb), is 
given by 



Qv t = 



\( w k \v b )\ 2 
<7i 



so that the total failure probability is 



Ipovm = V1Q1 + 



9i 



(5.4) 



(5.5) 



Choosing q\ to minimize the right-hand side, we obtain 
the expression for Qpovm given in Eq. H3.11JI . Again, we 
must have (2 < f]i < Ci f° r this expression to be valid. 

Now let us look at some choices for 771. If all of the 
functions are equally probable, then 771 = \/{Mb a i + 1), 
and the first von Neumann measurement is the optimal 
one. The measurement fails if the vector is \uik), but this 
event is so unlikely that it has a neglgible effect on the 
average failure probability. If 771 = 1/D, that is the prob- 
abilities of getting \wk) or a balanced- function vector are 
proportional to the dimensions of the subspaces in which 
they lie, then the POVM is optimal. If 771 = 1/2, so 
that we are equally likely to be given \wk) or a balanced- 
function vector, then the second von Neumann measure- 
ment is optimal. This illustrates how the best strategy 
is influenced by the a priori probabilities. 



VI. CONCLUSION 

If one is given a Boolean function that is promised to 
be either even or in W% , classically, in the worst case, one 
would have to evaluate it 2 n [(l/2) + (l/2 fe )] + 1 times to 
determine to which set it belongs. Using quantum infor- 
mation processing methods, one has a very good chance 
of determining this with only one function evaluation. 
This shows that Deutsch-Jozsa-type algorithms need not 
be limited to constant functions; certain kinds of biased 
functions can be discriminated as well. 

Unambiguous state discrimination is a procedure that 
is of fundamental interest in quantum information theory. 
Its only application so far has been to quantum cryptog- 
raphy. The results presented here suggest that related 
methods can also serve as a tool in the development of 
quantum algorithms. 
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Appendix 

We want to evaluate the sum in Eq. (|5.2(l . The vector 
\wk) is given by 

|iw fc ) = -^=(l,l,...l,-l,-l,...-l), (6.1) 



where the first D — (D/2 k ) places are l's and the fi- 
nal D/2 k places are —l's. Now consider a vector corre- 
sponding to a balanced function, which has to l's in its 
last D/2 k places, where < to < D/2 k , so that it has 
(D/2) - to l's in its first D - (D/2 k ) places. The overlap 
of \wk) with this vector is 



4m 



2 k ~ 1 D 



(6.2) 



The number of balanced functions of this type, C m , is 
given by 



r _( D/2 k \ ( D(2 k - l)/2 k 
m ~ I TO M {D/2) — TO 



so that the sum we have to evaluate is given by 
Sb Cn 



(6.3) 



ife-i 



D 



(6.4) 



where rj = (1 — Tji)/Mb a i is the a priori probability of 
each balanced function. 

In order to find Sb, we have to evaluate three types of 
sums. We shall discuss one is some detail, and simply 
give results for the other two. The first sum is 



D/2 k 

s = ^ G " 

m=0 



(6.5) 



This can be evaluated by noting that 

D/2" 



D/2 k = ^ ( D l 2 
TO 



1=0 



m— 

D(2 k - l)/2 k 
(D/2) - m 



c z (,6.6) 



Multiplying these two expressions together we find 

D/2 k D(l~2- k ) , , . 

(i+"=) D = e e r 

iv) — n l — n ^ / 



m=0 i=0 

D(l - 2- k ) 
(D/2) - m 



x' +m . (6.7) 



Comparing coefficients of x D / 2 on both sides of this equa- 
tion, we find that 



D 

■ s " 1 £>/2 
The remaining two sums are 

D ( D-l 



(6.8) 



D/2 k 

si = toC„ 



m=0 



2 k 1 P/2) - 1 



(6.9) 



and 



D/2 k 

S 2 = TO 2 C„ 

m=0 



2 k \2 k \(D/2)-2 



D_( D-l 
'2 k \(D/2)-l 



(6.10) 



Substitution of these expressions into Eq. I|6.4|l yields the 
result in Eq. I|5.3[> . 
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